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We prepare arbitrary patterns of neutral atoms in a one-dimensional (ID) optical lattice with 
single-site precision using microwave radiation in a magnetic field gradient. We give a detailed 
account of the current limitations and propose methods to overcome them. Our results have direct 
relevance for addressing of planes, strings or single atoms in higher dimensional optical lattices 
for quantum information processing or quantum simulations with standard methods in current 
experiments. Furthermore, our findings pave the way for arbitrary single qubit control with single 
site resolution. 
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I. INTRODUCTION 

Neutral atoms trapped in optical lattices form a 
promising paradigm of quantum simulation [1 and quan- 
tum information processing [2 . Numerous proposals in 
these fields require the ability to coherently address and 
manipulate atoms on a single site of the lattice. This 
ability, however, still poses a challenge in these systems. 
While detection of atoms in optical lattices with single 
site resolution has been reported [3 SI |5] O |7] , single site 
manipulation in lattices with a site separation in the opti- 
cal wavelength domain has so far only been demonstrated 
by removing atoms using a focussed electron beam [8 . 

A convenient method for spatially resolved coherent 
manipulation and detection of atoms using only global 
techniques originates from nuclear magnetic resonance 
(NMR). Originally developed for solid state systems, the 
basic tools and concepts have found their way also into 
quantum optics experiments using cooled trapped atoms 
[9 . Achievements include control of internal states on a 
micrometer scale [10] or robust global control of atomic 
samples [11]. In large systems, where a Bose-Einstein 
condensate is loaded into an optical lattice, periodic pat- 
terned loading was achieved by using superlattice poten- 
tials [T2j [13]. In an array of double wells, recently one 
well out of each double well system could be selectively 
manipulated by an optically induced effective magnetic 
field, while the qubit was stored in field insensitive inter- 
nal states [14]. 

Here, we demonstrate the use of NMR techniques to 
prepare arbitrary patterns of atoms in a ID optical lat- 
tice with single-site precision, which can be used as a 
starting point for quantum information processing, and 
can be extended to higher dimensional systems. These 
techniques overcome several restrictions posed on atom 
string preparation using moving optical lattices [TF , such 



* Electronic address: karski@uni-b onn.de| 
"^Electronic address: widera@uni-bonn.de 



as the limited resolution on the order of the beam diam- 
eter or the limited selectivity for closely spaced atoms. 
In principle, the NMR techniques can be realized with 
standard methods already used in current quantum gas 
experiments and can be extended to yield arbitrary quan- 
tum state manipulation of qubits with single site resolu- 
tion. We discuss the challenges of this extension in the 
context of our ID optical lattice. 

II. EXPERIMENTAL SETUP 
A. Cooling and Trapping of Atoms 

We capture neutral Caesium (Cs) atoms from the back- 
ground gas in a three-beam magneto-optical trap (MOT). 
The atoms are transferred into a far detuned standing 
wave optical dipole trap (ID optical lattice) formed by 
the interference of two counterpropagating laser beams 
with a wavelength of Aiat = 866 nm. By superposing 
the two traps for a duration of 500 ms and subsequently 
switching off the MOT, depending on the parameters, 
between a single and up to 100 atoms can be loaded into 
the optical lattice. For our parameters, atoms occupy a 
region of approximately 100 jiia length of the lattice (see 
Figj2|, corresponding to roughly 200 lattice sites. 

During the loading and fluorescence imaging phases 
(see below) the dipole trap has a depth of ks x 0.4 mK, 
whereas for manipulation of the atoms with microwaves 
in a magnetic field gradient it is adiabatically lowered 
within 50 ms to a depth of ks x 80 /iK. Here, in each 
dipole trap beam, approximately 20 mW of laser beam 
power are focussed down to a waist radius of 20 /im. 
This corresponds to vibrational frequencies of c^ax = 
27rx 115 kHz along the optical lattice axis (axial direction, 
2:-axis) and cjrad = 27r x 1.2 kHz perpendicular to the op- 
tical lattice axis (radial direction, x?/— plane). The atoms 
have a final temperature of /cb x 10 /iK, corresponding to 
mean vibrational quantum numbers of nax = 1-2 axial 
and nrad = 200 . 

The atomic sample is probed by fluorescence imag- 
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FIG. 1: (a) Zeeman splitting of magnetic sublevels in the 
6^5'i/2 ground state manifold of caesium. The degeneracy of 
the magnetic sublevels is lifted by an external guiding field due 
to the linear Zeeman effect. State |0) = |F = 4, mp = 4) and 
state |1) = |F = 3, = 3) define the states of the qubit. (b) 
Geometrical arrangement of laser beams and coils producing 
the magnetic fields needed for preparation and detection of 
arbitrary patterns of atoms. 



ing. We typically illuminate the trapped atoms between 
200 ms and 1 s by a near-resonant optical molasses us- 
ing the MOT laser beams, and image the fluorescence 
onto an EMCCD (electron multiplying CCD) camera. 
The particular exposure time is chosen depending on 
the required precision for the determination of atomic 
positions [6]. For our imaging system and the fluores- 
cence wavelength of 852 nm the diffraction limit is 1.8 /im 
^ 4 X Aiat/2. However, for a sparsely filled lattice, even 
in unresolved clusters of less than eight atoms, numerical 
post-analysis allows high-precision, real-time determina- 
tion of the atomic position with down to nearest neighbor 
distances [6]. 



B. Internal State Manipulation and Detection 

We initialize the atoms in state |0) = 
|F = 4,mi? = +4) by optical pumping with a de- 
polarized laser beam, which is resonant with the 
F = 4 ^ = 4 transition. Here F is the total atomic 
angular momentum and mp its projection onto the 
quantization axis, determined by a weak (|Bo| = 3G) 
guiding magnetic field along the z axis. An equally 
polarized repumping laser beam, frequency stabilized on 
the F = 3 ^ F' = 4 transition, is used to transfer the 
atoms back to the optical pumping cycle, whenever they 
decay to the F = 3 ground state. 



As information is stored in our experiment in the hy- 
perfine states |0) and |1) = \F = 3,mF = +3), we have 
to selectively probe the populations of these states. This 
is done using the so-called "push-out" technique [iF, 
which removes the atoms in F = 4 (including state |0)) 
from the optical lattice while leaving those in F = 3 unaf- 
fected. For this, we apply an intense laser beam operating 
on the F = 4 ^ F' = 5 transition, perpendicular to the 
optical lattice axis. The power Pp^gh ~ 
pulse duration r^^^^i — ^50 /is of the push-out beam are 
optimized so that its radiation pressure force overcomes 
the radial dipole force and quickly pushes the atoms in 
F = 4 out of the lattice within less than half a radial os- 
cillation period. Thereby, we largely prevent off-resonant 
excitations to = 4, from where the atoms can sponta- 
neously decay into the F = 3 ground state, which would 
cause erroneous detection of F = 3. We verified that the 
mean survival probability of atoms prepared in F = 4 is 
smaller than 1%, whereas for atoms prepared in F = 3, it 
is larger than 99%, imposing lower limits of the push-out 
efficiency. 

The two internal states |0) and |1) are coupled us- 
ing microwave radiation around 9.2 GHz. The microwave 
field is created by mixing a fixed frequency signal from 
a phase locked dielectric resonator oscillator (PLDRO) 
around 9.0 GHz with a tunable frequency from a sig- 
nal generator around 200 MHz in an upconverter. Both 
sources are locked onto a Rubidium clock. The result- 
ing sum frequency is amplified up to 12 W and directed 
onto the atoms using a commercial waveguide that sup- 
ports linear polarization, while the carrier and difference 
frequencies are suppressed by at least -30 dB. 

We thereby achieve Rabi frequencies of approximately 
1^0 = 27r X 60kHz for the |0) ^ |1) transition. The 
coherence time of the system has been measured by a 
Ramsey-type pulse sequence to be T2 ~ 200 /is and can 
be increased to T2* ~ 0.8 ms by a spin echo sequence. 



C. Magnetic Gradient Field 

In order to create a spatially varying resonance fre- 
quency along the optical lattice, we apply a quadrupole 
magnetic field Bquad(r, /). The field originates from two 
coils in anti-Helmholtz configuration, where the symme- 
try axis of the coils coincides with the optical lattice axis. 
At the center of the trapping region, the axial field gradi- 
ent along the z-axis is thus twice as large as the gradient 
in the x?/— plane. The magnetic field in the trapping re- 
gion including the offset field then reads 

B(r,/)= 1^ J j +B'{I) j , (1) 

where B'{I) denotes the magnitude of the gradient field 
along the 2;-axis for a chosen current / running through 
the coils. Due to the linear Zeeman effect, the |0) ^ |1) 
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transition frequency depends linearly on the modulus of 
B(r,/), yielding 



cJo(r,/) = AHFS+7|B(r,/)| 

= AHFS + ^0+^'(/)^+^P^ (2) 



where the gyromagentic ratio J = {3gs — 4:g4)/jjB /h ^ 
27r X 2.5 MHz/ G, = 7^o denotes the contribution of 
the guiding field, u;'{I) = jB'{I) the position-dependent 
shift and p = ^/x^^^y^ the radial distance from the 
axis of symmetry of the coils. Finally, |B(r, /)| has 
been approximated to second order in p, valid for 
{Bo + B\I)z)^ > B'{lfp^/A. Note that, due to mis- 
alignment, the optical lattice axis can be radially offset 
by a small amount po relative to the axis of symmetry. 
While the linear dependency of uJo{z,p,I) on the axial 
position z is still maintained, this offset imposes a signif- 
icant obstacle to single site addressability, as is shown in 
Section [mCl 

From a calibration measurement tracking the shift of 
the microwave frequency in the lattice when certain mag- 
netic gradient fields have been applied (see next section), 
we infer the position dependent frequency shift of 



27r/ 



(291 ±2) 



Hz 



(A/2) A 



(671 ± 3) 



Hz 
pmA 



(3) 



This frequency shift corresponds to a strength 
of magnetic field gradient along the lattice axis of 
5'(/)// = -(274± l)/iG/(/imA). For the maximum 
current of 45 A, we find a frequency separation for two 
neighboring lattice sites of uo' /{27r) x A/2 = 13 kHz, cor- 
responding to a magnetic field gradient of approximately 
120G/cm. 



III. SINGLE PARTICLE OPERATIONS WITH 
SPATIALLY VARYING RESONANCE 
CONDITION 

A. Position Dependent Spectroscopy 

When a magnetic field gradient is applied during the 
microwave pulse, the atomic resonance condition is ful- 
filled only in small regions along the lattice axis. In or- 
der to illustrate this, we apply a rather weak magnetic 
field gradient of B^ = 27.4 G/cm. For an atomic sam- 
ple which has been initialized to the |0) state, a subse- 
quent microwave pulse with rectangular amplitude enve- 
lope of duration tp^ise = 10 /is and a center frequency 
of 9.2 GHz is applied. As a consequence, the frequency 
spectrum with its characteristic sinc-type side-peaks is 
mapped onto the position distribution of atoms in |1). 
This is directly observable in the fluorescence images. 
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FIG. 2: (a) Superposition of 50 fluorescence images after a 
rectangular microwave pulse has been applied to a filled op- 
tical lattice in a magnetic gradient field to flip the internal 
state from |0) to |1) and a subsequent laser pulse has removed 
atoms in |1). (b) Vertically binned intensity distribution of 
the images shown in (a). The solid line shows a fit with the ex- 
pected sinc-like frequency profile. The insert shows the atom 
distribution without application of the microwave pulse and 
push-out laser. 



taken after the push-out has been applied, see Fig. [2j 
The averaged intensity distribution of several images thus 
directly illustrates the shape of the microwave pulse in 
frequency domain. Microwave spectroscopy in position 
space is therefore a useful tool, which immediately and 
quickly reveals information encoded in the frequency do- 
main with a high resolution, controlled by the strength 
of the field gradient. This method is much faster than 
a usual scan in frequency across the resonance in a ho- 
mogeneous field to map out the full spectrum. Further, 
using microwave pulses with a narrow-band spectrum, it 
allows us to monitor the evolution of the transition fre- 
quency arising from changes and drifts of experimental 
parameters from shot to shot. We stress that the shape of 
the microwave spectrum in position space is broadened 
due to the optical imaging when features in frequency 
space become smaller than the optical diffraction limit. 
A further broadening may exist due to a radial offset of 
the lattice axis with respect to the axis of symmetry of 
the coils. This effect, however, becomes only significant 
for narrow-band pulses with a high spectral selectivity as 
we discuss in Section BH C 41 



B. Patterned Atomic String Preparation 

The preparation of convenient initial atom configura- 
tions is an important first step for applications of neutral 
atoms for quantum information technology, such as the 
creation of entangled states through coherent collisions; 
the tailoring of atom strings for efficient interaction with 
the field of a high- finesse resonator; or the extraction of a 
selected plane or string of atoms from a Mott-insulating 
state of atoms |T7l [181 IE] • The concepts presented above 
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provide a toolbox which offers this capabihty: In presence 
of a magnetic field gradient, only those lattice sites re- 
main occupied after the state-selective push-out, at which 
atoms have been shifted resonant to the microwave tt- 
pulse. 

In order to generate a predefined pattern structure a 
pulse train of N successively applied 7r-pulses with dif- 
ferent frequencies coi is used instead of a single pulse. By 
incorporating the periodicity of the optical lattice using 



gN with z,j = l,...,iV,z^j, (4) 
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these frequencies define the pattern structure. To ensure 
a high selectivity in the position-dependant population 
transfer, we use the gradient providing the maximum 
frequency shift of cc;'(/)/(27r) x (A/2) = 13 kHz. Further- 
more, Gaussian 7r-pulses with a l/\/e spectral half-width 
down to cr^/27r = 6 kHz are employed rather than rect- 
angular pulses in order to suppress side lobes. For this 
maximum available frequency shift, these pulses are in 
principle capable of manipulating individual atoms with 
almost single site resolution. To increase the selectivity in 
the preparation of the pattern structures, we optionally 
use several iterations of state initialization, application 
of TT-pulses and the state-selective push-out. 

In Figure |3) two exemplarily generated patterns are 
shown: A string of eight equidistant atoms (Fig. [3]^a,b)), 
for instance, is ideally suited for a quantum register. Such 
a register has previously been implemented using atoms 
randomly distributed over the lattice sites ^ , in which 
case the position of each atom needed to be determined 
and fed back to the microwave source prior to the quan- 
tum state manipulation of the atoms. For a predefined 
string of atoms, this feedback is not required, since the 
frequency of each individual atom of the string is defined 
by the preparation sequence itself. 

Compared to rearranging atoms in crossed movable 
ID lattices [151 120] , our method allows us to prepare 
strings of atoms separated by an arbitrary number of 
lattice sites, even down to two (see Fig. [5|a,b) below). 
These small separations are essential for an efficient im- 
plementation of controlled collisions of two individual 
atoms ^21] [22] [23, ,24^ using, e.g., state-selective trans- 
port [251 [26]. 



Limitations 
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FIG. 3: Strings of atoms in an optical lattice: Greyscale im- 
ages show a single image acquired after application of the 
patterning sequence for (a) a predefined string of equidistant 
atoms with separations of 17A/2, and (c) for a predefined 
string of four "atom pairs" with an atom separation of 12 
lattice sites within each pair. Corresponding averages over 
50 acquired images are shown as false color images, (b) and 
(d), respectively. Atoms deviating from predefined positions 
are indicated by an arrow, the deviation by a boxed value. 
Missing atoms in (b) are mostly attributed to the limited 
preparation efficiency (see text). 



resulting in a filling factor of about 0.5, yielding a prob- 
ability to find an atom at one selected lattice site of 
Pa ^ 50%. This probability imposes an upper limit for 
the efficiency of generating the entire pattern structure: 
For a pattern of N atoms, the probability that all desired 
lattice sites are initially populated is given by pini = p^, 
yielding p\^\ ^ 0.4% for a string of 8 atoms. Since our 
atom detection is restricted to sparsely filled lattices, a 
detailed investigation of the loading process with high 
atom densities could not be performed so far. Further 
investigations aiming on the increase of the filling factor 
are required in the future. 



1. Filling factor 



2. Selectivity of microwave pulses 



The first restriction of our method is given by the ini- 
tial loading from a MOT. For high densities as in an 
optical lattice, light-induced collisions [27l [28] lead to 
high loss rates. In our case, any doubly occupied lattice 
site illuminated by the near-resonant molasses is quickly 
depleted |29]. The resulting steady-state distribution of 
atoms in the lattice will show only either one or no atom. 



For an atom at a predefined lattice site, the selective 
microwave operation, a 7r-pulse for instance, should fulfill 
two conditions. First, it should operate on the selected 
atom with high fidelity. The selectivity of this 7r-pulse, in 
turn, specified by its spectral width and its shape in the 
frequency domain, determines to what extent neighbor- 
ing atoms are affected by this pulse and whether neigh- 
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FIG. 4: (a) Sequence for preparation of pattern structures in- 
volving multiple application of the inner loop of state transfer 
and push-out (blue shaded blocks), (b) Microwave spectra in 
frequency domain for different numbers of applications M of 
the inner loop. Solid lines show Gaussian fits, from which 
(c) the dependency of the inferred maximum transfer effi- 
ciency -Pmax 

(M) (logarithmic plot) and (d) the 1/y/e spectral 
half- width auj{M) (double- logarithmic plot) on the number 
of applications is inferred. Both dependencies perfectly agree 
with the expected trend. 



boring atoms are completely removed from the lattice or 
remain trapped with a certain probability. For efficient 
preparation of pattern structures, this probability should 
by ideally zero. 

However, the selectivity cannot be infinitely improved 
by increasing the pulse duration. We observe experimen- 
tally that for pulse durations longer than 15 /is the max- 
imally achievable population transfer in the pulse center 
decreases together with the spectral width. We attribute 
this to decoherence stemming from, e. g. technical noise 
or inhomogeneous broadening due to radial motional dy- 
namics in combination with the differential AC Stark 
shift due to the lattice light [S^. For this reason, we 
employ a simple technique, which at least for the genera- 
tion of patterned structures presented in this work, effec- 
tively improves the selectivity of the Gaussian 7r-pulses. 
The underlying idea is to utilize the probabilistic and 
destructive nature of our state-selective push-out. 

Suppose that a single application of the microwave 
spectroscopy sequence with a Gaussian 7r-pulse provides 
a spectrum given by 



P\l){S) = Pmaxexp 



(5) 



with a certain maximum population transfer Pmax- 
Then, by repeatedly applying the inner loop comprising 
state initialization, application of the 7r-pulse and the 
push-out by a total of M times (see Fig. [4|, we expect 
again a Gaussian shaped spectrum 



Pmax GXp 



= (Pmax)'^exp 



M 



(6) 



2{(t^/VM) 



with a l/^/e spectral half-width ai^{M) = a^j^M, how- 
ever, with a reduced mgiximum population transfer 

Pmax(M) = (Pmax)"^ • 

In Figure [4j the scaling behavior of both quantities is 
exemplarily shown for a Gaussian pulse with gi = 20 /is, 
where the spectra have been taken in the frequency do- 
main rather than by position dependent spectroscopy. 
The scaling is perfectly reproduced by the measured data, 
for both quantities inferred from a fit. Comparing the 
performance of a sequence repeating the inner loop with 
a less selective 7r-pulse of = 15 /is by a total of two or 
three times with a single application of a more selective 
pulse of Gt = 20 /is, we conclude that, for the efficient 
generation of patterns, the repeated application is more 
advantageous. 

To investigate the efficiency, the resolution and pos- 
sible imperfections of our patterning method, we gener- 
ate a pattern of three pairs of next nearest neighbors, 
see Fig. [5) Each atom pair is separated by 16 lattice 
sites from another, to reliably determine the distances of 
the atoms using numerical post-processing as presented 
in Ref. [6 . The internal sequence core employs a pulse 
train of six subsequent applied Gaussian 7r-pulses with 
Gt = 20 /is. The inner loop is repeated by a total of 
M = 2 times. Consequently, we expect a maximum pop- 
ulation transfer of Pmax 

(2) = (71 ± 2)% at the predefined 
sites. The corresponding 1/^/e half-width in position 
space is cr2,exp(2) = (0.34 ± 0.01)A/2, specifying the se- 
lectivity region for a single atom. 

The entire patterning sequence is repeatedly recorded 
500 times. From the final fluorescence images of the 
atoms, each acquired with an exposure time of 800 ms, 
the positions and the distances of the atoms are deter- 
mined. For our parameters, the number of pairs cor- 
rectly prepared with the predefined distance is expected 
to be A/exp ~ 170 ± 17. In Figure l5] a histogram of the 
measured distances between atoms forming the pairs is 
shown. It reveals a Gaussian shaped distribution cen- 
tered at the predefined separation of two lattice sites. 
From the histogram, we infer A/'meas = 40 correctly pre- 
pared pairs. From the fitted width of the distribution 
c^dist we infer a selectivity region of a single atom as 

^meas = ^Tdist/ = (0.60 ± 0.06) A/2 . 

Both, the number of correctly prepared pairs and the 
1/v^ half- width of the selectivity region deviates from 
the expected values. These deviations can be attributed 
to the axial drift of the optical lattice and the radial offset 
of the lattice axis relative to the axis of symmetry of the 
coils, as discussed below. 
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FIG. 5: Pattern of three trapped atom pairs with nearest 
neighbor distance (a) in a single and (b) the average of several 
fluorescence images, similar to the patterns shown in Fig. [3] 
(c) Histogram of distances between two atoms forming the 
pairs. The solid line shows a Gaussian fit to the histogram, 
from which the 1 / half- width of the distances and the se- 
lectivity region of a single atom are inferred. 



S. Axial drift of the lattice 

The standing wave field of the optical lattice is suscep- 
tible to drifts and fluctuations of the optics, effectively 
changing the position of the potential minima over the 
course of a measurement. To infer the drift of the optical 
lattice relative to the imaging optics, we track the posi- 
tion of a single atom over a long time period from succes- 
sively acquired images with Is exposure time each, see 
Fig. [6^. We observe small fluctuations of the atom po- 
sition around an approximately linear trend, indicating 
a slow drift of the optical lattice of about lOnm/s. We 
attribute this drift to thermal expansion of the optical 
table and opto-mechanics of the lattice. 

Because the total measurement time is much longer 
than the time interval in which the lattice axially drifts 
over a distance of A/2, the effect of this drift on the spec- 
trum in position space can be described by a convolution 
equation 



I 4 

^|l),drift(^) = j jP\l){z - Z')dz' , 



(7) 



where we have included the fact that only drifts mod- 
ulo A/2 are distinguishable and relevant. By solving the 
convolution equation, we infer a "drift-free" l/>/e half- 
width of the selectivity region of (Jz^d^ = (0.52±0.07)A/2. 
This value explains part of the deviation between mea- 
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FIG. 6: (a) Drift of the optical lattice with respect to the 
imaging system over long times, (b) New coordinate system 
for an offset between the symmetry axes of optical lattice and 
gradient field, (c) Measured quadratic dependence of the po- 
sition of the microwave resonance frequency with increased 
offset ro along the z-direction. The offset is changed by ap- 
plying a homogeneous offset field along the ^direction. 



surement and expectation mentioned above. 



4' Radial offset of the lattice 

According to Eq. ([2| the transition frequency changes 
quadratically with increasing radial offset between the 
coil axis and the lattice axis. In the following, we calcu- 
late the effect of a radial offset of the lattice axis with 
respect to the axis of symmetry of the coils on the max- 
imum population transfer and the selectivity region of a 
Gaussian 7r-pulse as used above. 

To calculate the position dependent detuning, the ori- 
gin of the lattice coordinate system for r', which so far 
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was located at the center of the lattice site where the 
magnetic quadrupole field vanishes, is shifted by the off- 
set of the gradient field tq (see Fig. [gJd). For a position 
vector r in this new coordinate system, the spatial de- 
tuning is given by 

S{t\ To) = ujo{t' + ro, /) - cJo(ro, /) , (8) 

with a;o(r, /) from Eq. ([2|. 

In Figure [7) we show the calculated regions addressed 
by the pulse (addressed regions) for which |(5(r', ro)| < cTuj 
and |(5(r',ro)| < in presence of the field gradient 
for different values of radial offset pQ. Atoms located 
in these regions experience a minimum population trans- 
fer of Pmin > 0.6Pmax and Pmin > O.lSPmax, respectively, 
whenever the Gaussian 7r-pulse is applied. Because of the 
quadratical dependence of (5(r',ro) on the radial com- 
ponent, see Eqs. ([8| and the addressed regions are 
extending over neighboring lattice sites for increasing val- 
ues of pQ. This yields a decrease of the spacial selectivity 
of the Gaussian 7r-pulses. Furthermore, we observe a de- 
creasing intersection of the Gaussian thermal wave packet 
of a trapped atom and the addressed region. This leads 
to a decreased effective population transfer at desired lat- 
tice sites. We quantify both observations by calculating 
the expected position dependent spectra, i. e. the popu- 
lation transfer as function of axial displacement z' from 
the resonance position = for different radial offsets. 
For this, we average the population transfer specified by 
the spectrum of the Gaussian 7r-pulse (Eq. ([5|) over all 
position-dependent detunings, taking into account the 
axial and radial width of the atomic thermal wave packet 
in the trapping potential of the lattice. 

This yields again an approximately Gaussian shaped 
spectrum in position space, see Fig. [7j from which the 
effective maximum population transfer Pmax and the ef- 
fective l/>/e half- width in position space can be in- 
ferred. 

By comparing our measured spectrum to the calcu- 
lated effective spectra, we deduce a calculated offset of 
po = (64 ± 14) pm. For this offset, we expect a maxi- 
mum population transfer of Pmax = (56 ±8)%, and thus, 
an expected number of correctly prepared atom pairs of 
TVexp ^ 40 . This value agrees with the measured number 
of correctly prepared pairs stated above. 

Moreover, by changing the homogeneous offset field 
Bo, it is possible to infer and to reduce the offset between 
magnetic gradient and optical lattice axes, see Fig. 
For this we record the spatial position of the microwave 



resonance along the gradient as a function of the radial 
offset pq. This offset is increased by applying a homo- 
geneous field along the x'-direction. We observe the ex- 
pected quadratic dependence of the resonance position. 
Adjusting the magnetic field offset to the measured value 
of the parabola maximum, the offset along the scanning 
direction can be removed. Repeating this for the orthog- 
onal direction we are able to remove this offset completely 
from our system. 

This effect is also important for the selection of 2D 
planes from a perfect 3D Mott-insulting state with unity 
filling. We estimate that for a spherical symmetric cloud 
of 25 pm diameter, for ideal alignment and with our lat- 
tice, gradient and pulse parameters, a 2D plane of 2500 
atoms could be prepared by applying the inner loop of 
Fig. [4] twice. Here, approximately 3% of the atoms stem 
from neighboring planes. 



IV. CONCLUSION AND OUTLOOK 

We have created arbitrary patterns of atoms trapped 
in a ID optical lattice with single-site resolution using 
magnetic resonance techniques. Such patterns can be 
used as initial configurations for various purposes, in- 
cluding collisional interaction between nearest neighbors. 
The mechanisms limiting the resolution have been iden- 
tified as being technical in nature and can be overcome. 
Specifically, the optical lattice potential can be actively 
stabilized to avoid long term drifts, the magnetic field 
gradient can be increased and the microwave pulses can 
be shaped by methods of optimal control to increase the 
selectivity of the pulses while preserving maximal trans- 
fer probability. In a further extension of this method, 
arbitrary single qubit operations with single-site resolu- 
tion - even in massively parallel operation - may become 
possible by taking into account off-resonant phase shifts 
on other qubits. 
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